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Section 1.0 
INTRODUCTION 


Propagation delay through the Ionosphere Is one of the error sources 
for which a correction must be made for future satellite altimeters, with 
satellite to mean sea surface height accuracies required at the few centi- 
meter level. The maximum correction amounts to approximately 15 cm. at 
14 GHz. Due to the large variations between daytime and nighttime Ionospheric 
electron densities, a correction would be required even for precision pro- 
filing missions for which absolute height accuracy Is not required. 

Several approaches are possible, Including measuring the crrrectlon 
through the use of a dual frequency altimeter, and the use of various types 
of Ionospheric modeling. This report considers the two modeling techniques 
which are considered most promising for application to near real-time 
altimetry processing, with a recomnendatlon made as to the technique giving 
superior accuracy on a global scale. 
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Section 2.0 
IONOSPHERIC MODELS 

Two types of Ionospheric models have been proposed as a basis 
for making height corrections for satellite altimeter data. The basic 
characteristics of these two approaches are summarized below. 

2.1 NASA Model 

The "NASA" model Is an empirical worldwide model of the Ionosphere, 
and exists In several different versions. What will be referred to here 
as the NASA model [1] Is presently r irt of the operational Goddard 
Trajectory and Determination Syste.u (GTDS) at Goddard Space Flight Center. 

An earlier version Is also In operation at the Naval Surface Weapons Center 
at Dahlgren, Virginia. 

The NASA model assumes a profile shape as shown In Figure 1. This 

profile Includes a number of key parameters, with the maxImiAn electron 

density and height of the maximum density among the most critical. Typical 

global maps of height of maximum density (H^^) and critical frequency 

inox 

(corresponding to maximum electron density) are shown In Figures 2 and 3. 

All parameters have latitude, longitude and diurnal dependence, which 
are driven by dally solar flux (f^g y) values (plus averaged values over the 
previous year) and a set of M3000 coefficients available from NOAA/Boulder. 

A typical global map of total electron content, to which the Ionospheric 
propagation correction Is proportional, Is shown In Figure 4. 

With the Input of dally flux values, this technique gives corrections 
which are accurate to within approximately 20-25% [2]. Since the largest 
lonospheMc propagation correction (at 14 GHz) Is ~15 cm. for a maximum 
solar activity period, the NASA model errors could reach 4 cm. For less than 
peak solar activity, the error would be considerably smaller and the RSS 
error would probably be on the order of 1-2 cm. even for maximum solar 
activity, since nearly half of an average pass would be In nighttime for 
which the correction error Is near zero. 
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Figure 1. The Exponential, Parabolic and Bl-Parabollc 
Ionospheric Profile Used as a Basis for the 
NASA Ionospheric Model 
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Figure 2. World Map of Height of Maximum Ionospheric Electron Density 
Hmax at 12 h U.T. on January 15, 1969. Contours Every 20 km. 
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Figure 3. World Map of Ionospheric Critical Frequency fOF2 in MHz 
at 12 h U.T. on January 15, 1969. Contours are Plotted 
Every 1.0 MHz. 
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Figure 4. World Map of Integrated Vertical Electron Content N as 
Predicted by NASA Ionospheric Model In Units of 10* 
el/sq. meter at 12 h U.T. on January 1969. Contours 
Every 5 x 10^‘ el/sq.m. 


2.2 JPL Faraday Rotation Model 


The Faraday Rotation or "dPL" model has been used at the Jet Propulsion 
Laboratory, Pasadena, California. It Is based on measured Faraday rotation 
measurements from one or more ground stations to a satellite. These measure* 
ments can then be converted to Integrated election density measurements along 
the slant range from the ground station to the satellite (outside the Ionosphere). 
These measurements are then extrapolated [3] to other parts of the world as 
functions of time, spacecraft geomagnetic latitude, and local solar zenith 
angle. In the Implementation of this model for Seasat altimeter data cor- 
rections [3], two ground stations were used, Armidale, Australia for the 
southeren hemisphere and Goldstone, California for the northern hemisphere 
(both stations using the ATS 1 satellite). 

Based on various comparisons made during Seasat data evaluation. Including 
comparisons with corrections computed from the dual -frequency Geocelver data, 
the accuracy of the Faraday Rotation method has been estimated [3] to be 
3 cm. The corrections are, of course, very good near the measured Ionospheric 
path but can have errors on the order of 50% at geographical locations far 
from the measurements. 
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Section 3.0 
MODEL RECOMMENDATION 


Of the models discussed In the previous section, the NASA solar flux 
model* Is recomended for the following reasons: 

1. Since the Ionosphere Is a highly variable and Inhomogeneous medium, 
it does not extrapolate well beyond about 800 km. [4]. Thus, a limited 

set of Faraday rotation measurements Is not sufficient to provide cor- 
rections on a global basis with as high degree of accuracy as can be 
obtained from the global empirical model. 

2. Should something happen to the Faraday rotation data flow, the 
extrapolation technique then degrades rapidly. The global model, even 
without solar flux updates, degrades much more slowly. 


* The author would like to acknowledge that he Is hardly an unbiased 
observer with regard to Ionospheric models, having been deeply Involved In 
the development of what Is referred to here as the NASA Ionospheric model. 
It Is believed, however, that the arguments and evidence In support of the 
empirical modeling approach are overwhelming. 
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Section 4.0 

IMPLEMENTATION AND CHECKOUT OF NASA MODEL 

The NASA model In the STDS progrem already contains the NOAA M3000 
coefficient set. These coefficients consist of 36 sets of coefficients, 
each covering a 10 day time span, and have had seasonal and solar flux 
dependency removed. They are thus permanent coefficients and need not be 
periodically replaced. In operatlo.i, the data set Is stored on disk, with 
only a limited set read Into core as needed. 

The model Is currently operated as a subroutine and Its Implementation 
as a part of another system would require the removal of some linkage 
with 6T0S. But the subroutine package Is largely self-contained and the 
required modifications should be relatively minor. 

Model maintenance requires that dally solar flux values be continually 
provided. As previously Indicated, a set of values for the previous year Is 
required. 

Although the NASA model presently makes corrections for non- vert leal 
paths and for angle as well as ranging measurements, there would be minimal 
Improvement In efficiency or reduction In core requirements by making the 
program correct only altimeter data. It Is thus reconmended that the 
program be used basically as Is (and as a black box). Extensive checkout 
of a transferred program could be made through comparisons of Intermediate 
printout of the various critical parameters, as well as the computed 
corrections, and comparison with computed quantities on the base computer. 
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